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(54) Method for quantification of an analyte 

(57) The invention provides a new method for the 
quantification of an analyte, wherein in a first step, the 
analyte is amplified by an amplifying agent. Subse- 
quently, the amount of amplification product is deter- 
mined as a function of time and the extremae of first, 
second or n th order derivatives are calculated. The 
obtained values can then be used to calculate the initial 
concentration of the analyte. 
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Summary of the Invention 

[0009] The present invention describes methods for quantitatively analyzing either an amplif iable or a self replicat- 
ing system in such a way that the amount of amplification or self replication product is measured continuously. In 

5 accordance with one embodiment of the present invention a method for quantifying an analyte is provided. The method 
comprises the steps of contacting the analyte with an amplifying agent, and amplifying at least one predetermined locus 
of the analyte. The amount of amplification product is then determined as a function of reaction time and the first, sec- 
ond or n th order derivative of said function is calculated wherein n is a natural number. The maximum, the zero value 
or the minimum of said derivative is then determined and the initial concentration of the analyte is calculated from said 

10 maximum, zero value or minimum. The disclosed method can also be used to quantitatively analyzing the growth rate 
of organisms and thus measure the effect of a modified environmental condition on the growth of an organism. 



Brief Description of the Drawings 
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Fig. la: 
Fig. lb: 
Fig. 2a: 
20 Fig. 2b; 

Figs. 3a-3j: 
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Fig. 4: 

Fig. 5: 
30 Fig. 6: 



Measurement of fluorescence signal with FRET-HybProbes. 

Linear regression curve calculated with the threshold method according to Example 1 . 
Measurement of fluorescence signal with FRET-HybProbes. 

Linear regression curve calculated with the second derivative maximum method according to Example 
1. 

Plots showing fractional cycle numbers as calculated in six fold replicate by determining the maximum 
first derivative, minimum second derivative, and maximum second derivative versus log 10 of template 
copy number in a SYBR Green experiment according to Example 2. Savitzky Golay filter vertical win- 
dow parameters were set as follows Fig 3a. 2/3/2; Fig. 3b, 2/4/2; Fig. 3c, 3/3/3; Fig. 3d. 3/4/3; Fig 3e. 
4/3/2; Fig. 3f. 4/4/2; Fig. 3g. 6/3/2; Fig. 3h. 6/4/2; Fig. 31. 6/5/2; Fig 3j. 6/6/2. 

Graphic representation of the fractional cycle number determination by quadratic estimation of the 2nd 
derivative maximum. 

Graphic representation of a logistic growth curve representative of the growth of a population of cells. 
Graphic representation of the growth curves generated for a population of bacteria ceils cultured rep- 
resentative of the growth of a population of cells cultured in the presence of substance A (A), substance 
B (B) substance C (C) and the control (T). 



Detailed Description of the Invention 

35 

[0011] The invention comprises methods for quantitatively analyzing either an amplif iable or a self replicating sys- 
tem in such a way that the amount of amplification or self replication product is measured continuously Continuous 
measurement In a PGR reaction, for example, may also mean that measurements are performed once each cycle. The 
method of the present invention is especially applicable for embodiments, wherein during one phase of the amplification 
40 reaction or self replication (usually starting from the beginning or alternatively after an initial lag phase) the amount of 
the target product increases progressively Moreover, the method Is applicable for ennbodiments wherein after said 
exponential phase, the rate of amplification decreases. 

[0012] By using appropriate mathematical algorithms, the generated data can then be used to set up a mathemat- 
ical function, which discloses the smoothened kinetics of signal generation at different time points. In order to obtain a 

45 function, from which derivatives can be calculated, it is possible for example to perform a polynomal fit thereby creating 
a data set which reflects the smoothened kinetics of the originally measured raw data, it has been proven to be advan- 
tageous for the present method, if the polynomal fit Is calculated by a Savitzky Golay filter (Press et al., in: Numerical 
Recipes inC, the art of scientific Computing, Cambridge University Press, second edition, p. 650 ff). From such a func- 
tion, first, second or even higher orders of derivatives can be calculated and the maximal, minimal or zero values of said 

£0 functions can be determined. The obtained data are characteristic parameters for each amplification or self replication 
and can be used either to quantify the original amount of amplified entity or the influence of compounds or environmen- 
tal factors (such as temperature, pH. electromagnetic radiation, or nutrient composition) which are present during 
amplrf ication or self replication of the sample. 

[001 3] The majority of amplif iable or self replicating systems work in a way that amplification or self replication over 
55 a certain period of time is obtained exponentially, which is usually the case at the very beginning of the reaction. How- 
ever, in case of growth curves of microorganisms or cell populations, for example, there may exist a lag phase at the 
beginning, which at a later time point may enter an exponential phase. However, this does not affect the applicability of 
the invention for such embodiments. 
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^,uences, which are present in ^^^'^""f ^ ^.^^^^^^ add analytes is the polymerase 

aDPropriate amplilication methods. The most common of basic applications is given ,n S,e- 

S nLction (PGR), which is well-known for persons sMl«nnt^^^^^^^ New Jersey, p. 55-79 (1998)). 

bS . Molecula! Diagnosis of infectious ^'seases (ed^ Re s^^^^^ ^^^^.^ ^,,,3 e.^^ple NASBA 

Within the scope of the invention, however, ^^.^^^dl J^^^^^^^^ °NA. single-stranded DNA or any 

(MaleK et al.. WO 91 102814). T»>e ^rf t nuc.e.c aad r^^^^^^^^ ^^^^^^ be performed, 

type of RNA. in the latter case, prior to PGR. a cDNA syntnes « described, wherein the analyte 

rS.15] in accordance with one embodiment. « ^^"^^^^ an amplifying agent (and all com- 

s a Slynucteic acW. The method comprises ^^^/JJt^f theTnS a^^^^ at leas?a part of the analyte. Typ- 
pounds Which are necessary to amplify at least a pan "J ^ nucTeic acid sequence only the sequence 

fc^ly only a limited region of the analyte is amplrf.ed ^^^^Jl'^;.^!^^^ then be detemiined as a function 

by the two primers will be be calculated, wherein n is a r^ura. 

:ir;rrc;:=^^^ 

reaction. . . rietermination of amount of amplification product as a function of Ijm^ 

[0016] in the context of this invention, the deterrnination ot^a n,easured. Instead, it is sufficerrt to 

Lsiot mean, that the absolute amour, o a ^^^^^ 

determine relative signals, since ^^'ative determination al eaay^ concentration, 
to the invention may then be used to ^^t^^'^^ ^^//^^.J^j^^^;^^^^^ is obtained by a polymerase chajn reac- 

1001 7] in one embodiment of the invention the f the present invention, the amplication 

ion ?nd the analyte to be quantified is ^^''^^^''ZJ^^^^^^^ - a nucleic acid, the amplification product 
product is typically detected by means of ^'^o-^escence^ \Amen the a^^^ y ^^g^^ binding fluorescent 

Cbe^etTcted b^ a double-stranded DNA bin^ng^^^^^^ 

dye Alternatively, the amplrfication P"""^"^^ "^^^ .f j!^^^ 

p^es. which have sequences « J ^° ^^^^^^^^ 

range of at least 1 0 <=°"«9"°"f J^^"^^, '^^^^^^^ are hybridized to the amplification product. Fluores- 

each labeled with a fluoresced entty ^^^ ^In't^k^ o.ace between the two fluorescent entitles. _ 
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analyte and said binding entity IS formed 

b) determining the amount tunction. wherein n is a natural number. 

organism o. • cell poBuWoo. oompnamo 

„ g,o«», s.id mic«.rg=nl.m » "I' l»!»^;^^^^fp!;^T at~«on ol 9.<»* «™ 

t, 3a„™9 tha "» numb=, 01 ''^ ™J°SS^S2KSor. «h.,«n n is a nataml ^.J^ 

rrlC'Sr.r:;:,^?:=rr**a.».s.^eo,«..eo,o.sa.»™.^^ 

^^^^c..^.«.so..»..o...™«a,e,*as^...ecn»«„o.,™o^-a™o™ 
S** can .N, b, pro.aad in »n,i,«l '^"^^^^^^^^'^^Bl^l*^- wherein ,» amoTtotiOO 
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PGR In this case, the reaction may be limited by the amount of DNA polymerase, the concentration of amplification 
primers or the concentration of nucleoside triphosphates. Alternatively, progressive amplification may also be limited by 
product inhibition; e.g. in the case of PGR by reannealing of single stranded PGR products to form double stranded 
products that cannot hybridize to the primers. 
5 [0022] According to the new invention, a data driven algorithm is used to select data for further analysis. In order to 
obtain an appropriately smoothed kinetics, a mathematic filter algorithm can be used. In case of nucleic acid amplifica- 
tion, preferably, a Savitzky Golay filter may be applied (Numerical Recipes inG, the art of scientific Gomputing, Gam- 
bridge University Press, second edition, p. 650 ff). 

[0023] Upon applying a Savitzky Golay filter, each measured kinetic point is replaced by a smoothened point, cal- 
10 culated from a window including points to the left and to the right of the neighborhood. From these points, a polynomial 
of n th order is calculated for every measured value; with the exception of those values determined at the very beginning 
and the very end of the reaction, where due to lack of points to the left or to the nght, such a calculation is impossible. 
Accordingly, vertical window parameters are defined as a/n/b, wherein a is the number of points taken from the left, n is 
the order of the polynomial and b is the number of points taken from the right. 
15 [0024] From the smoothened kinetic values, the maximum, minimum or zero value of the first, second or n th deriv- 
ative is calculated. Determination of such extremae provides suitable set points for the definition of an unique and reli- 
able fractional cycle number characteristic for each kinetic, which reflects the initial concentration of the analyte. 
[0025] It has been proven to be particular advantageous to use a Savitzky Golay filter having the vertical window 
parameters 3/3/4 in combination with the calculated maximum of the first derivative, or, alternatively, a Savitzky Golay 
20 filter with the vertical window parameters 6/2/4 in combination with the calculated maximum of the second derivative in 
order to obtain the most reliable data. 

[0026] For each smoothened kinetic point, derivatives are calculated by methods known in the art. For example the 
first derivative for each smoothened value y for a given cycle number x of a polymerase chain reaction maybe calculated 
as f'(x) = y(x)-y(x- 1) The second derivative can then be calculated from the first derivative in the same manner Alter- 
25 natively, at a given cycle x the polynomial fit from the Savitzky Golay filter can be easily evaluated for first and second 
derivatives of y and x. 

[0027] Subsequently, in order to determine the extremae from said derivatives, a further polynomal fit may be 
applied, for which a number of three to five points are taken into account, which are located around the presumed 
extremum. For example Fig. 4 shows the fractional cycle number determination of the second derivative maximum 
30 using a quadratic fit to five points. The second derivative values were determined by a sliding window polynomial fit (a 
Savitzky Golay filter). 

[0028] The possibility of measuring the kinetics of an amplification reaction, which is a prerequisite for the new 
invention, has become enormously facilitated since there are instruments and methods available, wherein the genera- 
tion of the amplification product can be measured continuously by spectroscopic detection principles, for example by 
35 means of fluorescence. An example for a suitable instrument is described in detail in Wittwer et al., Biotechniques 22, 
No. 1. 176-181 (1997). 

[0029] Several detection formats based on target dependent fluorescent signaling have been disclosed, which ena- 
ble continuous monitoring of the generation of amplification products (reviewed in Wittwer etal., Biotechniques. Vol. 22. 
No, 1, 130-138, 1997). These detection formats include but are not limited to: 

40 

1. Use of fluorescent double-stranded DNA recognizing compounds 

Since the amount of double stranded amplification product usually exceeds the amount of nucleic acid origi- 
nally present in the sample to be analyzed, double-stranded DNA specific dyes may be used, which upon excitation 
with an appropriate wavelength show enhanced fluorescence only if they are bound to double-stranded DNA. Pref- 
45 erably. only those dyes may be used which like SYBR Green I, for example, do not affect the efficiency of the PGR 

reaction. 

2. Increased fluorescence resonance energy transfer upon hybridization 

For this detection format, two oligonucleotide hybridization probes each labeled with a fluorescent moiety are 
50 used which are capable of hybridizing to adjacent but non overlapping regions of one strand of the amplification 

product. Preferably, one oligonucleotide is labeled at the 5' end and the second oligonucleotide is labeled at the 3' 
end When hybridized to the target DNA. the two fluorescent labels are brought into close contact, such that fluo- 
rescence resonance energy transfer between the two fluorescent moieties can take place. As a consequence, the 
hybridization can be monitored through excitation of the donor moiety and subsequent measurement of fluores- 
55 cence emission of the second acceptor moiety. 

In a similar embodiment only one fluorescently labeled probe is used, which together with one appropriately 
labeled primer may also sen/e as a specific FRET pair (Bernard et al.. Analytical Biochemistry 235, p. 101-107 
(1998)). 
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3. Tag Man principle 

In order to detect the amplification product, a single-stranded hybridization probe is used, which is labeled with 
a fluorescent entity, the fluorescence emission of which is quenched by a second label on the same probe which 
may act as a quenching compound During the annealing step of the PGR reaction, the probe hybridizes to its tar- 
get sequence, and. subsequently, during the extension of the primer, the DNA polymerase having a 5'-3'-exonucle- 
ase activity digests the hybridization probe into smaller pieces, such that the fluorescent entity is separated from 
the quencher compound. After appropriate excitation, fluorescence emission can be monitored as an indicator of 
accumulating amplification product. 

4. Molecular beacons 

Similar to the Taq Man Probes, a molecular beacon oligonucleotide is labeled with a fluorescent compound and 
a quencher compound, which due to the secondary structure of the molecule are in dose vicinity to each other. 
Upon binding to the target DNA, the intramolecular hydrogen bonding is broken, and the fluorescent corrpound 
located at one end of the probe is separated from the quencher compound, which is located at the opposite end of 
the probe (Lizardi et al., US 5, 11 8,801). 

[0030] Usually, the hybridization probes as disclosed have sequences which are completely identical with or exactly 
complementary to the sequence of the analyte. However, it is also within the scope of the invention, if the probes con- 
tain one or several mismatches, as long as they are capable of hybridizing to the analyte under appropriate hybridiza- 
tion conditions. In any case, it has been proven to be particular advantageous, if the sequence identity or 
complementarity is 100% over a range of at least 10 contiguous residues. It has also been proven to be advantageous, 
if the length of the probe does not exceed 100 nucleotides, preferably not more than 40 nucleotides. 
[0031] In case of quantification of an analyte according to the principle of using external standards, a calibration 
cun/e has to be generated. When the analyte to be quantified is a nucleic acid which is amplified by PGR. known 
amounts of the analyte are amplified and the intensity of fluorescent signal is determined as a function of cyde number. 
After smoothening of the kinetics by a mathematical fit. one or more extremae of the first, second or n th derivative are 
calculated. This enables a correlation between the original target concentration and the fractional cycle number of a 
determined extremum. Subsequently, determination of unknown analyte concentrations may be performed. 

Example 1 

Amplification of a CyclophillinA sequence by PGR and subsequent detection with HybProbes: 
Comparison with a previously disclosed algorithm 

[0032] From a linearized plasmid carrying the GyclophillinA sequence, was subjected to serial dilutions to produce 
concentrations of 10^. 10^ 10^, andlO^ copies per 2 jul were prepared. In order to obtain reliable statistical data, these 
dilutions were prepared in 6 fold replicates. In a total volume of 20 |il. for each sample a PGR reaction was set up: 



1 x LCDNA Master mix hybridization probes (Roche Mol biochem. Catalogue No. 2015102) 
4 mM MgCl2 

1 fiM GycA1 fonward primer (pos. 52-69) 

1 ^iM GycA2 reverse primer (pos. 471 -493) 

200 nM 3* Fluorescein CycA hybridization probe (pos. 420-439) 

200 nM 5' LC RED640 CycA hybridization probe (pos. 444-470) 

10^. 10^. 10^, and 10^ copies of CycA plasmid. 



[0033] Nucleotide positions of primers and probes as indicated corespond to nucleotide positions as given in 
Haendler et al. (Embo Journal 6. p. 947-950 (1987)). 

[0034] According to manufacturer s protocol, the samples were filled in 1 .5 mm capillaries (Roche Mol. Biochem. 
Cat. No. 1909 339). Subsequently, the capillaries were inserted into a Roche Diagnostics GmbH LightCycler LC 32 and 
subjected to a PGR reaction according to the following thermocycling protocol: 
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Temp. (LC) 


Time (sec) 


Rampe rate ( JC/s) 


Acqu. 


Cycles 


Denaturation 


95 


120 


20,0 


None 


X 1 


Denaturation 


95 


0 


20.0 


None 


X45 


Annealing 


55 


10 


10.0 


Single 


X45 


Elongation 


72 


15 


3,0 


None 


X 45 



[0035] After each cycle, according to manufacturer's protocol, fluorescence was analyzed for each sample in two 
different channels, wherein channel 1 corresponded to a bypass of 520 +/- 20 nm and channel 2 corresponded to a 
15 bypass of 645 +/- 20 nm. Subsequently, the ratio of the values obtained from channel 2 versus channel 1 was calculated 
in order to normalize volume differences due to pipetting errors and exclude alignment and detection errors. 
[0036] Raw data obtained in Example 1 were processed in two different ways: 

A) According to the previously disclosed method, a manually chosen threshold level was set after appropriate back- 
20 ground subtraction of raw fluorescent data. For each sample, the intersection between the threshold line and the 

function of fluorescent signal versus time was calculated with the following algorithm: Based on a regression cal- 
culation, two measured points above a previously defined noise band were used to define a log linear phase by 
means of linear regression. Within the log linear area, a threshold line was set manually, and the intersection 
between the threshold line and the regression lines corresponding to each tested copy number were calculated. In 

25 order to gain information on the validity of the method, the calculated crossing points were used to create a plot 

indicating cycle number versus calculated log 10 of initial concentrations of the target DNA (see Figs, la and Ib). 
From this plot, a linear regression was determined, which enabled the calculation of mean values, standard devia- 
tions and variation coefficients as indicators for the accuracy and precision of the method. Results are shown in 
Table 1, wherein concentrations and crossing points were calculated with the fluorescence threshold method 

30 according to Example 1 , and the abbreviations are as follows: mean values (Mittelwert) standard deviations 
(STDWN) and variation coefficients (CV(%)). 

B) According to the new invention, the same raw data were smoothened by a polynomal fit using a Savitzky Golay 
Filter having the vertical window parameters 6/2/4 without previous background subtraction. From the obtained 

35 function, the maximum of the second derivative was calculated and used as an indicator for the initial template con- 
centration. Again, from the obtained data, a plot indicating the fractional cycle number versus initial concentrations 
of target DNA was created (see Figs. 2a and 2b). Mean values, standard deviations and variation coefficients were 
calculated as described above, results of which are shown in Table 2 wherein concentrations and crossing points 
were calculated with the second derivative maximum method according to Example 1 , and the abbreviations are 

40 as follows: mean values (Mittelwert) standard deviations (STDWN) and variation coefficients (CV(%)). 

[0037] As can be seen by comparison of results from Tables 1 and 2, the statistic values either obtained by method 
A) according to prior art, or obtained by method B) according to the invention yielded about the same degree of accu- 
racy and precision. However, since the method according to the new invention is independent from the actual amount 
45 of generated fluorescent signal and data processing is not user driven, it provides a superior advantage over the previ- 
ously disclosed quantification methods. 

Example 2 

50 Amplification of a TFN-alpha sequence by PGR and subsequent detection with SYBR Green: 

Comparison of calculating initial concentrations from the first derivative maximum, the second derivative maximum and 
the second derivative minimum. 

[0038] The experiment was carried out according to Example 1 with the following modifications: 
55 [0039] A plasmid carrying a TNF-alpha genomic sequence (Shirai et al.. Nature 313. p 803. 1985) was subjected 
to serial dilutions to produce sample concentrations of 10^. 10^, lO'^, 10^ and 10^ copies per 2 lJ. 6 fold replicates were 
prepared for each sample concentration. 

[0040] In a total volume of 20 nl. each PGR reaction contained: 
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1 X LCDNIA Master mix SYBR GreenI (Roche Mol. brochem. Catalogue No. 2015099) 
3 mM MgCl2 ' 

1 TNF-alpha forward primer (Roche Mol. Biochem. Cat. No. 1 989 626) 
1 nM TMF-alpha reverse primer (Roche Mol. Biochem. Cat. No. 1 989 626) 
10^, 10^. 10^. lO^and 10^ copies o1 TNF-alpha cDNA containing plasmid. 



[0041] The applied thermocycling protocol was identical to Example 1. After each cycle of amplification. SYBR 
Green fluorescence was measured according to the manufacturers' protocol in LightCycler channel 1 . corresponding to 
a bypass of 520 +/- 20 nm. 

[0042] According to the new invention, the raw data were smoothened by different polynomal fits using a Savitzky 
Golay Filter having the vertical window parameters 2/3/2. 2/4/2, 3/3/3, 3/4/3. 4/3/2. 4/4/2. 6/3/2. 6/4/2, 6/5/2. and 6/6/2. 
From the obtained kinetics, the maximae of first derivatives, second derivatives and the minimae of second derivatives 
were calculated for all six replicates of a given analyte concentration and used as an indicator for the initial template 
concentration. 

[0043] Based on these data, cycle numbers were plotted against the loglO of different copy numbers of template 
DNA. Subsequently, a regression line was calculated. As can be deduced from Fig. 3a-3j. the results of these calibra- 
tion performance plots demonstrate the linear relation between the calculated cycle number and the log 10 of analyte 
concentration, proving that principally determination of different extremae by using different types of vertical window 
parameters may be used for determining initial analyte concentrations. 

[0044] In order to determine preferable embodiments regarding the selection of appropriate Savitzky-Golay filters 
with respect to the first and second order maximae. the data from the different plots were compared mainly by taking 
the following 3 criteria Into account: 

a) variation coefficient CV (%) as determined for the 6 replicates representing 100 copies of the analyte 

b) slope of the straight regression line, as an indicator for the separation potential 

c) separation potential for discriminations between 10^ and 10^ copies of analyte. calculated by the robust statisti- 
cal Wilcoxon Test wherein p values are determined, for which a number of less than 0.05 is regarded as significant. 

[0045] The results of this analysis are summarized in Table 3. 



Derivative 


RIter-Parameter 


Slope 


Separation i^otential 
IO0IOO copies 


CV% at 100 copies 


1st 


2/3/2 


-3.29 


.0042 


1 


1st 


2/4/2 


-3,33 


.0038 


1 


1st 


3/3/3 


-3.29 


.0040 


0.85 


1st 


3/4/3 


-3.29 


.0041 


0.85 


1st 


4/3/2 


-3.30 


.0040 


0.76 


1st 


4/4/2 


-3.29 


.0043 


0.93 


1st 


6/3/2 


-3.29 


.0041 


0.85 


1st 


6/4/2 


-3.30 


.0037 


0.81 


1st 


6/5/2 


-3.29 


.0039 


0.9 


1st 


6/6/2 


-3.31 


.0041 


0.91 


2nd 


2/3/2 


-3.42 


.0042 


1.26 


2nd 


2/4/2 


-3.27 


.0036 


7.69 
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(continued) 



Derivative 


Filter-Parameter 


Slope 


Separation Potential 
IO0IOO copies 


CV% at 1 GO copies 


2nd 


3/3/3 


-3.37 


.0035 


0.67 




o/4/o 


-0.4o 


.UUoo 


1 


2nd 


4/3/2 


-3.40 


.0037 


1.18 


2nd 


4/4/2 


-3.41 


.0037 


1.78 


2nd 


6/3/2 


-3.39 


.0043 


0.96 


2nd 


6/4/2 


-3.39 


.0044 


1.09 


2nd 


6/5/2 


-3.40 


.0040 


1.09 


2nd 


6/6/2 


-3.32 


.0041 


1.69 



[0046] Regarding separation potentials, the data shows that using different filters prior to calculation of the first 
derivative maximum, as well as prior to calculation of the second derivative maximum, there exist similar separation 
potentials with respect to criterium b) and sufficient separation potentials for alt filters tested with respect to criterium c). 
[0047] For determining the maximum of the first derivative, comparatively low CV(%) values have been obtained for 
all tested filter parameters. For determining the maximum of the second derivative, vertical window parameters of 6/2/b 
are identified as preferable, wherein b Is either 3, 4, or 5. 

Example 3 

[0048] Logistic growth curves can be used to model population growth (Fig. 5). Such curves have a beginning lag 
phase, an exponential growth phase, and a plateau phase. Bacterial growth follows this kind of a curve. The effect of a 
stimulating or inhibiting compound on growth rate can be determined by following the relative number of organisms over 
time. If the curve shifts to the left, the compound is an inhibitor, if it shifts to the right, the compound stimulates growth. 
The potency of inhibition or stimulation is related to how far the curve shifts. 

[0049] For example, a growth media (LB broth) in a spectrophotometer cuvette is inoculated with a bacteria {E. 
coli), stirred and kept at 37DC. Every 5 minutes, bacterial growth is monitored by the light scatter through the cuvette 
at 500 nm. As the number of bacteria increase, the absorbance increases, giving a relative measure of the number of 
bacteria. The second derivative maximum is calculated to give a fractional cycle number of the control sample (T). 
[0050] In parallel reactions, substances A, B, and C are added to separate cuvettes with media and inoculated as 
before. Their grovi^h is monitored over time as above. Curves such as shown in Fig. 6 will be obtained. As indicated in 
the demonstrative data of Rg. 6, substance A is a more potent inhibitor than substance B, while substance C is a stim- 
ulator. The potency of inhibition can be quantified by the magnitude and direction of the shift in fractional cycle number: 

C J- C B = magnitude of inhibition of substance B 

C -p- C ^ = magnitude of inhibition of substance A 

C -p- C Q = magnitude of inhibition of substance C 

[0051] In addition to measuring the turbidity of the media (as measured at 500 nm) to assess growth, it should be 
clear that many other monitors of growth can be used, including pH shifts of the media or consumption of a nutrient. 
The only requirement is that a relative measure of population size is obtained sufficient to calculate derivatives. 

Claims 

1 . A method for quantification of an analyte. comprising 

a) contacting said analyte with an amplifying agent: 

b) amplifying at least one predetermined locus of the analyte. 

c) determining the amount of amplification product as a function of reaction time, 

d) calculating the first, second or n th order derivative of said function . wherein n is a natural number, and 
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e) determining the maximum, zero value or minimum of said derivative 

0 emulating Lm sa,d max.mum, zero value or minimum the inrtial concentra^on of the analyte. 

2. The method of Cairn 1 , wherein during one phase of the amplification reaction the amount of ampmcation product 
5 increases progressively. 

3. The method of claim 2, wherein after said progressive phase, the rate of arriplHication decreases. 

4. The method of claim 3, wherein said derivative is calculated by a mathematical fit. 

5. The method of cteim 1 . wherein the amplification product is detected by means of fluorescence. 

6. The method of claim 4, wherein the arrplif ication product is detected by means of fluorescence. 
,5 7. The method of claim 1 . wherein the analyte is a nucleic acid. 

8. The method of claim 6, wherein the analyte is a nucleic acid. 

9. The method of claim 8. wherein the amplification is obtained by a polymerase chain reaction. 
1 0 The method of claim 4. wherein the amp.«ication is obta.ned by a polymerase chain reaction and the amplication 

product is detected by a double-stranded DNA binding entity. 

tiguous residues, 
the two fluorescent entities. 

13. A method for analyzing the effect of a compound on the growth rate of a microorganism or a cel. population, com- 
prising. 

pound. 
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14. A method for quantification of an analyte. comprising 



b) contacting said analyte with a signal generating binding entity under conditions where a complex between 
said analyte and said binding entity is formed 
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FIG. 3A 
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FIG. 3B 
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FIG 3C 
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FIG. 3E 



EP1 041 158 A2 



45- 
40- 

35- 

30 

CYCLE NUMBER 

25-1 

20 



15- 



10 ^ 



MAX_1DER 

MIN 2DER 




MAX 2DER 



n 1 1 1 1 r 

2 3 4 5 6 7 

LOG10(CONC) - COPY NUMBER 



17 



FIG. 3F 
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FIG. 3G 
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FIG. 3H 
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FIG. 3J 
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